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We investigate the nonlocal conductivity calculated from the time-dependent Ginzburg-Landau 
equation. When fluctuations of the vector potential are negligible, the high-temperature (Gaussian) 
and low-temperature (flux-flow) forms of the uniform conductivity within an ab plane {(Jyy) are 
essentially identical. We find to what extent the nonlocal conductivity shares this feature. The 
results suggest that for pure samples in these regimes the length scales of the nonlocal resistivity, 
Pyy{y — y',z — z'), remain short-ranged in the y and z directions in contrast to the assumptions made 
by the hydrodynamic modelling of the multiterminal transport measurements. On the other hand, 
the resistivity is seen to have a long length scale in the x direction pyy(x — x'). The implications for 
the interpretation of recent experiments are discussed. 
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The conductivity is nonlocal if the current at a site r 
is determined not only by E(r) the electric field at r but 
also by the fields at other sites E(r'). It is always nonlo- 
cal on some microscopic scale, but the term "nonlocal" is 
generally reserved for cases in which the length scale in- 
volved is much longer than some underlying microscopic 
one. A nonlocal conductivity implies a nonlocal resistiv- 
ity (though the length scales on which each varies may 
differ). Since the resistivity of type II superconductors 
will be influenced by the motion of flux lines, it seems 
reasonable to expect nonlocal effects. In recent exper- 
imjants on heavily twinned samples of YBCO, Safar et 
aZ.EJ have observed the onset and enhancement of nonlo- 
cal conductivity as the temperature is lowered. 

With experiments on clean samples, samples with 
columnar defects and films underway, it would be useful 
to sort out the basic form of the conductivity. But the 
vortex-liquid state of a type II superconductor is a prob- 
lem with many length scales and identifying which is/are 
relevant to a particular measurement is by no means triv- 
ial. 

In the Ohmic (linear response) regime, nonlocality im- 
plies that the (dc) conductivity (T^y(r,r') ^ a^„S{r — r') 
or in momentum space: cr^,y(k) ^ const. Husg_and 
Majumdam explored a "hydrodynamic" approachH'El, in 
which the small-k expansion of the conductivity is trun- 
cated as follows: 



crpt.(k) = 0-^^(0) 
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In the example they work out in detail to explain fea- 
tures of the multiterminal transport measurements, they 
include only one S: Syzzy It models the "tilt viscos- 
ity" — which measures the influence of pancake vortices 
moving in one ah plane on those mojang in another when 
they tr|ayEl with different velocities .□ The hydrodynamic 
picturg3'Q began by viewing vortices as the analogues of 
polymers in a melt, concepts like viscosities naturally fol- 
lowed. Thus this picture would seem most appropriate 



in the London limit where the core size and magnetic 
screening length A are shorter than the distance between 
vortices, i.e. where vortices are most polymer-like. It 
is unclear, however, what fraction, if any, of the vortex- 
liquid regime is described by such a hydrodynamic ap- 
proach. 

Stability requires (T^^(k) to be a positive definite ma- 
trix; for the hydrodynamic form of the conductivity, it 
implies for instance that Syzzy and Syyyy be positive. 
When these conditions are met, one finds that the re- 
sistivity decays with a length scale of order S/a), 
which is in turn the length|.acalc found in the electpic-field 
and current distributions .□ Recent calculationsQ-B start- 
ing from the time-dependent Ginzburg-Landau (TDGL) 
equation find that many of the S"s do not have the req- 
uisite sign — at least for clean samples at high tem- 
peratures. We study the nonlocal conductivity within 
the ab plane ayy(k) in the lowest Landau level (LLL) 
limit, which is valid not close to Hc2{T) as is some- 
times stated in the literaturg but at lower temperatures 
in the vortex-liquid regimeJZl When fluctuations of the 
vector potential are negligible, the uniform conductivity 
Gyy (0) is known to interpolate smoothly between its high- 
temperature (Gajissian) form and its low-temperature 
(fiux-flow) form.Q We extend these results to the non- 
local conductivity and argue that Syzzy and Syyyy do 
not change sign as the temperature is lowered, suggest- 
ing that there may be a substantial region where the 
hydrodynamic approach does not apply. 

Our starting point is the TDGL equation for the su- 
perconducting order parameter 5* in a constant magnetic 
field along the z axis in the Landau gauge: 



ie* Bx^ 



vl/(r,i)+6|vl/(r,i)p*(r,t) = 77(r,i), (2) 
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where the {ri(r,t)) = and has (5-function correlations 
{rj*{r, t) r,{r\ t')) = ^ 5[v - r') 5{t - t'). (3) 

At high temperatures, the nonhnear term (oc 6) in (|^) 
can be ignored. The "states" of the hnear operator col- 
lapse into Landau levels 



(4) 



where are the harmonic-oscillator eigenstates and 

£ = {h/e*By/^ the magnetic length. The corresponding 
energies are 



En{kz) ~ h^kl/2'mc + a + nhujQ, 



where loq = e* B/mab is the cyclotron frequency and a — 
a + the energy of the LLL. 

Some length scales are already apparent; clearly (. is 
an a6-plane length scale related to the distance between 
flux lines. Two energy scales, a and ?iwoj appear in 
En{kz)', from these one can construct two c-axis length 
scales: — jlm^aY^'^ and 4 = {^Im^ujof-I"^ = 
{'niabl'>ncY^'^L The former is the mean-field c-axis cor- 
relation length; the latter a c-axis analogue to the mag- 
netic length. The lengths, I and ^c^ depend only on the 
magnetic field; while grows as the temperature (T) is 
lowered. 

The Kubo formula can be used to determine the con- 
tribution of the superconducting fluctuations to the dc 
conductivity (the Aslamazov-Larkin term) : 
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The superconducting currents are given by: 
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Each involves a a and a (gauge-invariant) deriva- 
tive, so when performing the average in (|6|), one "con- 
tracts" *(r) with ^'*(r') and likewise **(r) with *(r') 
resulting in two correlators. In diagrammatic form the 
result is a loop, the first term in Fig. 1(a). 
The Gaussian result for o-yy (0) is 
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We will not reproduce the full expression for the Gaussian 

ajfv(k-) here as they appear in Ref. [5]. Terms in the 
small-k expansion are similar to that above in Eq. (|^). 
When a <C cjo, these sums are dominated by their first 
terms, yielding: 



(k)- 
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eiki + o{k'') 
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where $o =h/2e and (|^o|^), the equilibrium average of 
the fluctuations of the order parameter in the LLL, is 



ksT m. 
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First note that Syyyy and Syzzy are negative (non- 
hydrodynamic) and that the length scales multiplying 
(5) ky and k^ are magnetic lengths. On the contrary, Syxxy 



is positive and the associated length scale ^± oc uj^lja 
grows as T is lowered. So the first candidate for a growing 
length scale in the resistivity Pj,y(r — r') lies within the 
ah plane, a somewhat unexpected result as the picture 
of moving flux lines is more suggestive of a long c-axis 
length scale. The phase coherence length and the length 
over which density-density fluctuations decay have pre- 
viously been-|identified as T-dependent length scales in 
the ah planeO; recent Monte-Carlo simulations find that 
they grow in the same way.Ej It will prove interesting to 
follow the renormalization of to determine whether it 
is distinct from these or not. 

Next note that this high T result for cry?/(0) is 
very reminiscent of the flux- flow conductudty ct// = 
ri/<^QB (where is a viscosity coefficient). til The flux- 
flow conductivity is derived-from translating an unpinned 
Abrikosov flux-line latticolj, and its validity depends on 
whether one can neglect thermal fluctuations in the vcc- 
tor potential. Its similarity to the expression for (jyy 
in eq. (|^) occurs since lowering T has minimal effect 

on ay^y{Q). This feature is specific to the ah conductiv- 
ity; it is not shared by the c-axis conductivity ai^z{0). 
It does extend, however, to the nonlocal conductivity 
(Ty^y {kx = Q,ky,kz), though it no longer holds if kx 7^ 0. 
We exploit this latter point to argue that Syyyy and Syzzy 
do not change sign as T is lowered and that the associated 
length scales remain short. 

So what is this property? In the calculation of ay'^y (0) 
the momentum operators occurring in the definition of 
the currents (|^) act like creation (or annihilation) op- 
erators on the VP's and 4'*'s. They raise (or lower) the 
"states" with the following important outcome: of the 
two correlators that result, one is in the n*'' Landau level 
the other is in the {n + 1)*''. In the a <C limit, the 
contribution with one n — Q correlator and one n = \ 
correlator dominates. (See Fig. 1(a).) 
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FIG. 1. (a) Diagrams which contribute to the ab plane con- 
ductivity in the LLL limit, (b) Example of a diagram which 
is "down." (c) Diagram with vector-potential fluctuations 
which have not been included. 

Next recall that the eigenfunction expansion of the cor- 
relator C(r, r'; t') = (^{r, t + t')^*{r', t)} is 

C(r,r';t')oc/ / ^ WC(rO ^,r.„(..)|.| 

From here we see that in the LLL limit the time integral 
in the Kubo formula becomes 



dt' e-r[-Bo(p)+-Bi(p-fc)]t' 



(12) 



where p is the z component of momentum running 
through the n = channel and p — k that through the 
n = 1 channel. (See Fig. 1(a).) Note that Eq is much 



smaller than Ei and we can drop it from the integral ( 12 ) . 
This last step is equivalent to using the static (equilib- 
rium) n — correlator and yields the factor of (IV'oP) in 



The main point is that one and only one of the corre- 
lators is in the n — level and that it can be replaced 
by a static n = correlator. We now extend these ar- 
guments to the nonlocal conductivity. Consider Syzzy 
which featured in the hydrodynamic modellingcl of the 
multiterminal transport measurements. The only modi- 
fication in calculating Syzzy is that it requires two addi- 
tional derivatives with respect to z. As these derivatives 
do not affect the Landau-level structure, the above argu- 
ments carry through. In fact, we take further advantage 
of the disparity in "masses" (a <C wq) by insisting that 
the external momentum kz be sent through the n — 1 
channel. (This choice does not affect the outcome, just 
makes it more apparent.) The internal momentum inte- 
gral is dominated by poles associated with mass a, and 
the terms involving kz and wq are left essentially intact. 
The full kz dependence in the LLL limit is: 



[1 + klel/2] 



(13) 



The structure is simply that of [Ei{kz)]~'^ with a — > 0. 

More detailed calculations reveal that Syyyy has a sim- 
ilar property. This time the extra derivatives (with re- 
spect to y) have some effect on the Landau-level struc- 
ture: while one correlator is in the n — level, the other 
is in cither the n = 1 or n = 2 level. The important 
point is that there is still no contribution with two n ~ 
correlators. This last point applies to all terms in the ky 
expansion; the series resums to give: 



'yy V"-!/ 



{ky 
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In fact, all the terms in the small-k expansion of 
(Tyy{ky, kz) (so loug as kx — 0) share this property, and 
the combined ky and kz dependence is: 



'^yy {ky,kz) 



(n + l)(fcg£V2)" 
^,n\[{n + l) + k^Jl/2Y 



On the other hand, the extra derivatives in Sy 



(15) 



(with 



respect to x) have a more crucial effect, producing a 
nonzero contribution when both correlators are in the 
n = state. The outcome, as seen previously, is that 
Syxxy is positive (hydrodynamic) and the long length 
scale arises. 

We want to go beyond the Gaussian result, i.e. con- 
sider the effect of the nonlinear term in the TDGL equa- 
tion. One question is: To what extent is the simple loop 
affected? In the region of the phase diagram dominated 
by the LLL, the answer is "not much," as already hinted 
at by the similarity of the Gaussian and flux- flow results. 
Another question is: Where is this region? The Gaus- 
sian calculation suggests that the LLL approximation is 
valid near the Hc2{T) line (a = 0), but renormalization 
effects will change this! A very thorough discussion of 
the region of validity of the LLL appears in Refs. [7]. 

Our first attack on the problem of including the non- 
linear terms, i.e. renormalization effects, will be the 
Hartree-Fock (HF) approximation. The theory remains 
effectively Gaussian but a d = a + 26 (| -001), giving the 
following self-consistent equation for a = d + huJo/2: 
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b ksT m, 

2l/27r £27j 5,1/2- 
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Now we ask not where a is small but where a is small to 
establish where the LLL approximation is valid. The 
renormalized mass a grows small as the bare mass 
a grows large and negative, at temperatures below 
Hc2{T). Of coursg, renormalization effects beyond the 
HF approximationQ ought to be considered, but concern- 
ing the issue of where the LLL approximation is valid, 
they are only refinements; the essential features are cap- 
tured by HF. 
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To go beyond the HF approximation to a'ifj {k^ = 
0,ky,kz), we expand perturbatively around it. Contribu- 
tions to the TO*'' order require 2m additional contractions 
(to correlators and to propagators). The pivotal issue is 
in what Landau levels do these new correlators and prop- 
agators lie. li a ^ ujQ, then any contribution in which 
they are in higher Landau levels would be "down" (have 
a significantly lower weight) compared to those in which 
they are all in the LLL. Any insertion that disrupts the 
n > 1 correlator (as in Fig. 1(b)) automatically results 
in an extra higher Landau-level contraction. The only 
way in which all new terms are in the LLL is if they 
only dress up the n = correlator. This result is shown 
schematically in Fig. 1(a). 

As a consequence, in the LLL approximation, <Jyy{kx = 
0,ky,kz) still involves only two correlators: a fully renor- 
malized n = correlator C*'"' (r, r', i') and an n > 1 
correlator untouched beyond the HF approximation. In 
the LLL limit (7^°^ (r, r', i') has exactly the same depen- 
dence on X, x' , y and y' as its "bare" version; thus, the 
/cy-dependence of dyy is exactly the same as that in ex- 
pression ([l^). Moreover, a repetition of the arguments 
following eq. (|l2|) shows that the n — Q correlator can 
still be considered static, that the external momentum k^ 
can still be sent through the n = 1 channel and that the 
internal momentum integral is still dominated by poles 
associated with mass a. So as before the terms involving 
the external momentum kz and mass remain essen- 
tially intact. The upshot is that the /c^-dependence of 
dy^y is exactly the same as that in expression ( [Tsl ) and in 
particular Syzzy and Syyyy remain negative. Only (IV^oP) 
gets renormalized. 

The perturbative arguments above (summarized in 
Fig. 1(a)) do not hold when vector-potential fluctua- 
tions are no longer negligible. Then we would include 
vertices involving a 5", a and an A^, where is the 
fluctuating part of the vector potential. These vertices 
lead to diagrams like that shown in Fig. 1(c), where the 
wavy line represents an A-A correlator. In the Landau 
gauge, vertices associated with Ay come with a factor of 
loq which is large in the LLL limit. They also involve a 
(gauge-invariant) derivative and so the ^ and ^E"* at such 
vertices are in different Landau levels (see Fig. 1(c)). 
Diagrams with A-A correlators disrupting the n = 1 ^'- 
^ correlator (as in Fig. 1(c)) are not down compared to 
those with A-A correlators dressing up the the n = cor- 
relator. Hence the previous arguments no longer apply if 
vector-potential fluctuations are important. 

The effects of disorder will be considered in a future 
publication. Briefly, point disorder enhances the conduc- 
tivity in the ah plane; however, the c-axis length scale 
appears to remain short, ^c■ (In a nutshell, the external 
momentum kz can still be sent entirely through an n = 1 
channel in the leading diagrams.) On the other hand, 
in a case with correlated disorder, such as columnar de- 
fects or twin planes, the c-axis length scale grows as T is 
lowered, suggesting an increased likelihood for nonlocal 



effects. 

In summary the central result is that in the LLL ap- 
proximation, which holds overj-a considerable portion of 
the H — T diagram for YBCO}!\, the nonlocal conductiv- 



ity cTyy [ky, kz) has the form given in Eq. (15) with spe- 
cial cases given by Eqs. (13) and ( p| ) with only (|i/'oP) 
having any temperature dependence. In particular, the 
length scales multiplying ky and k^ in these expressions 
are short and temperature- independent; in addition, the 
signs of Syyyy aud Syzzy rcmalu negative (nonhydrody- 
namic) throughout the regime. Hence one would expect 
in experiments on clean samples with geometries that 
probe chiefly ayy to see essentially local behavior, as may 
be borne out by recent experiments.113 
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